Conical Flow induced by Quenched QCD Jets 



J.Casalderrey-Solana, E.V. Shuryak and D.Teaney 

| Department of Physics and Astronomy 

Q ! State University of New York, Stony Brook, NY 11794-3800 

O 



(N 
3 



oo 



o 



X 

5-H 



E-mail: casalderOtonic .physics . sunysb . edu, shuryak@tonic .physics . sunysb . edu, 
dteaney@tonic .physics . sunysb . edu 



Abstract. Quenching is a recently discovered phenomenon in which QCD jets created in heavy 
ion collisions deposit a large fraction or even all their energy and momentum into the produced 
' matter . At RHIC and higher energies, where that matter is a strongly coupled Quark-Gluon 

Plasma (sQGP) with very small viscosity, we suggest that this energy /momentum propagate as 
■ a collective excitation or "conical flow". Similar hydrodynamical phenomena are well known, 

^ ' e.g. the so called sonic booms from supersonic planes. We solve the linearized relativistic 

, hydrodynamic equations to detail the flow picture. We argue that for RHIC collisions the 

direction of this flow should make a cone at a specific large angle with the jet, of about 70°, 
and thus lead to peaks in particle correlations at the angle A(j> — n ± 1.2 rad relative to the 
large-pt trigger. This angle happens to match perfectly the position of the maximum in the 
angular distribution of secondaries associated with the trigger recently seen by the STAR and 
PHENIX collaborations. We also discuss briefly possible alternative explanations and suggest 
some further tests to clarify the mechanism. 
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a: 

_C • 1. Introduction 

Jet quenching is an important phenomenon, predicted in a number of papers and recently 
observed at RHIC 0. The main research has been so far related to the calculation of the energy 
loss dE/dx of the fastest parton, and focused on a kind of tomography of the produced matter. 
In this work, however, we focus on a different question: Where does the energy of the quenched 
jets go? We point out that as the hydrodynamical description of sQGP excitations works well, 
it should be used to predict how collective flow would develop, after a local deposition of energy 
and momentum. As we will see below, it will be a coherent source of sound waves in conical 
form. 

In what follows, we would like to treat two different types of energy losses separately : 
(i) the radiative losses, producing relativistic gluons in the forward direction; and (ii) the 
scattering /ionization losses, which deposit energy and momentum directly into the medium, 
as well as radiative losses of gluons at rather large angles (see Discussion section). Such gluons 
are rather soft and are promptly absorbed by the medium. 

The study of radiative losses were started in [3], then corrected for the destructive 
interferences (the so-called Landau-Pomeranchuck-Migdal or LPM effect) in i.. For a recent 
brief summary, see also 0. Although it is the dominant mechanism for quenching hard {z ~ 1) 
fragmenting hadrons, in our problem the second type of losses are more important because the 
primary parton and the radiated gluons all move with speed close to the speed of light and are 
treated as one object. 



Figure 1. A schematic picture of flow created by a jet going through the fireball. The trigger 
jet is going to the right from the origination point (the black circle at point B) from which sound 
waves start propagating as spherical waves (the dashed circle). The companion quenched jet is 
moving to the left, heating the matter and thus creating a cylinder of additional matter (shaded 
area). The head of the jet is a "nonhydro dynamical core" of the QCD gluonic shower, formed 
by the original hard parton (black dot A). The solid arrow shows a direction of flow normal 
to the shock cone at the angle 9m, the dashed arrows show the direction of the flow after the 
shocks hit the edge of the fireball. 



Elastic energy losses were first studied by BjorkenpP, while those due to "ionization" of 
bound states in sQGP were recently considered by Shuryak and Zahed (BJ. These mechanism 
deposit additional energy, momentum and entropy into the matter. (Like for delta electrons 
in ordinary matter, this excitation kicks particles mostly orthogonal to the jet direction.) It is 
their combined magnitude, dE/dx = 2 — 3GeV/ fm, the one we will use below. Even at such 
loss rates, a jet passing through the diameter of the fireball, created in central Au-Au collisions, 
may deposit up to 20-30 GeV, enough to absorb the jets of interest at RHIC. 

Let us start our discussion of associated collective effects by recalling the energy scales 
involved. While the total CM energy in a Au-Au collision at RHIC is very large (about 40 
TeV) compared to the energy of a jet (typically 5-20 GeV), the jet energy is transverse. The 
total transverse energy of all secondaries per one unit of rapidity is dE±/dy ~ 600 GeV. Most 
of it is thermal, with only about 100 GeV being related to collective motion. Furthermore, 
the so called elliptic flow is a ~ 1/10 asymmetry and therefore it carries energy ~ lOGe^ 
which is comparable to that lost by jets. Since elliptic flow was observed and studied in detail, 
we conclude that conical flow should be observable as well. (In order to separate the two, it 
is beneficial to focus first on the most central collisions, where the elliptic flow is as small as 
possible.) 

Fig^explains a view of the process in a plane transverse to the beam. Two oppositely moving 
jets originate from the hard collision point B. Due to strong quenching, the survival of the trigger 
jet biases it to be produced close to the surface and to move outward. This forces its companion 
to move inward through matter and to be maximally quenched. The energy deposition starts 
at point B, thus a spherical sound wave appears (the dashed circle in FigQ)- Further energy 
deposition is along the jet line, and is propagating at the speed of light, till the leading parton 
is found at point A at the moment of the snapshot. 

As is well known, the interference of perturbations from a supersonically moving body (such 
as a supersonic jet plane or a meteorite) creates a conical flow behind the shock waves. Similar 
flow was discussed in RefsJT] for shocks in cold nuclear matter, in which compression up to QGP 



production takes place. Unfortunately, experiments have shown that nuclear matter is too dilute 
and dissipative to make shocks: but discoveries made at RHIC let us be more optimistic in the 
sQGP case we consider now jH]. 

The angle 9m defined in the figure is given by a simple geometrical condition: the distance 
traveled by the jet during the proper time interval r is AB = ct, while the one traveled by the 
wave is 

CB = [ c s (t)dt =>- cosd M = — ( c s dt . (1) 

J0 CT Jo 

We use the speed of sound since, apart from the region close to the head of the object [S], the 
shock waves are weak and thus they move with the speed of sound (we use units in which c=l.) 

The region near the head of the jet, which we will refer to as a "non-hydrodynamical core" , will 
not be discussed in this work. Let us just mention that near point B, where the jet is produced, 
it consists of an "undressed" hard parton only. However, as it is constantly emitting gluons, 
which emit new ones etc., the whole shower is a complicated nonlinear phenomenon which should 
obviously be treated via the tools of quantum field theory [TOj. As found in the multiplicity 
of this shower grows nonlinearly with time, so eventually the core may become a macroscopic 
body, providing a large perturbation of the matter. Prom the hydrodynamical point of view, its 
size is limited from below by the dissipative "sound attenuation length" T s = (4/3)ry/(e + p), 
with T) being the shear viscosity. 

A shaded region in Fig^ consists of "new matter" related to the entropy produced in the 
process dS/dx, which can be calculated only with dissipative dynamics in the near zone, which 
we do not attempt in this work. Matter is expected to get equilibrated soon, and thus the radius 
of the cylinder R c of new matter can be estimated as dS/dx = s(T)irR^.. As is well known, a 
constant size cylinder does not emit any sound. 

2. Linearized Hydrodynamics 

The hydrodynamical equations we use are simpler than those used in Refs 7_ since we will 
consider that the total energy-momentum density deposited by the jets is a small perturbation 
compared to the total energy of the medium. This allows us to linearize the problem. The 
approximation breaks down close to the jet, where we will not describe the hydrodynamic fields. 
We will use cylindrical coordinates with x and r parallel and perpendicular to the jet axis, 
respectively. We will also assume that the perturbed medium is homogeneous and at rest. 

In the linearized approximation we define the following quantities in terms of perturbations 
of the stress energy tensor: 

e = ST 00 , g* = 5T 0i . (2) 
The remaining non-zero components of the stress tensor are: 

TV = <*eSV 1— (ay) , (&9>) = &9 i + dfg i -ls i ia i g i . (3) 

e +p \ ' \ I 3 

Here we have used 2| = c\ and recall that in the linearized approximation the velocity field is 

v l = eo+ p i e o and po are the energy density and pressure of the unperturbed medium. 

The energy and momentum conservation equations d^5T^ u , can be written in Fourier space, 

where we observe that if we define g = g-^x + ffr (where L ant T stand for longitudinal and 
transverse respectively), the linearized hydro equation decouples as: 



d t e + ikg L = , 
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k 2 gT = . 



(5) 



The system of equations (jU describes sound waves (propagating modes). Equation © is the 
diffusion equation and is not propagating. Only the sound waves will form the Mach cone. 

3. Initial conditions 

The initial conditions are set by the process of thermalization of the energy and momentum lost 
by the jet. This complex process should take place at distances of order T s from the production 
point. As T s is also the minimal size of the liquid cells, we will simply consider that there is a 
variation of T^ v at the position of the particle. 

If we consider an infinitesimal displacement dto of the high energy particle (moving with the 
speed of light) in the x direction at the point (to,to,0), we need to specify the infinitesimal 
variation of the fields edt (t = to,x), gdt {t = to,x). Due to the symmetries of the problem, the 
most general expression for those functions is: 



As argued before eo, go and gi are some functions with characteristic scale a ~ r s . The exact 
expression for these functions depends on the details of the thermalization process. As we do 
not know these details, we will need to make certain assumptions about the different functions 
that appear in (JHJ). In general, one may consider two different scenarios of matter excitation: 

Scenario 1 in which local deposition of energy and momentum is described by the first terms 
in (6), eo and go- This imposes a constraint in jHJ): 



where 4— and 4- are an input. In this case g\ = 0. These simple initial conditions, however, 
excite the "diffuson mode" and thus be discarded later. 

Scenario 2 in which the excitation is due to the gradient term g\ in (6), with eo = and 
go = 0. The deposition of energy and momentum in this case comes from the second order 
effects. The normalization of such solutions are better done at large distances, via the energy 
flow and the momentum flow through a large cylinder. As we show below, this generates the 
conical solution with sound excitation. 

The empirical fact that the second scenario describes the data and the first one does not 
obviously provides some insights about the excitation mechanisms, but we would not speculate 
about them now. 

4. Solutions 

In order to solve Q and @ with the initial condition @, we will calculate the corresponding 
kernels, that is, we will set the functions appearing in © to be <5-functions. Following the 
standard procedure, we convolute these kernels with some parametrization for eo, go, and g\\ 








T^(x) = J d 3 y (K£ (x - y) e (y) + K% (x - y) g (y) + (x - y) 9l 



(y)) 



(8) 



and integrate over the trajectory of the particle. 

Scenario 1. In this case we only consider point like excitation in eo, go- Keeping terms to 
order 0{T 2 S ) we find the following 4-vector kernels: 



i=± 



(9) 
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Where r = \x — x(to)\ is the distance from the observation point to the jet position and we have 
defined the functions: 

1 (r±c s t) 2 3 1 £— 

P ± (r s ) = -_e-^i^ t P( T S )= e . (11) 

V2vrr s t 2 ^/27r|r s i 

The second term in (|10j) is the (Fourier transform of the) so called diffusion mode (or diffuson) 
(JHJ). It creates a dissipating flow of matter comoving with the jet that remembers the location 
of deposition of momentum. As we will show below, this type of excitation leads to a forward 
peak in the final spectrum filling up the Mach cone (related to © and the first term in (|1U[|). 

These kernels should be convoluted with the initial distributions. We take a simple 
parametrization (Gaussian) that fulfills the requirement Q: 

. dE e 2^ dpx e —2 

eo( * = *°' s) = ^°(W)372 > 3o(t = t , X) = ^10-^^8 . (12) 

Scenario 2 Following the same procedure we find the 4-kernel: 

K = E (^%W.),#^arPi<r.)) , (i3) 



i=± 



where the definitions are the same as in @ and (|lUj) . Note that as claimed, we now find a 
propagating solution, the sound wave. 

For example, let us assume a simple parametrization for the initial source: 

gi(t = t , x) ~ e 272— _ ( 14 ) 

In Fig|2]we show the velocity field along the jet axis generated by the high energy particles 
according to both models of the initial conditions. In Fig. Eta) we can observe the effect of the 
diffuson as a maximum in the velocity along the x, axis, where the jet propagates. In Fig. Efb) 
such maximum does not exist, as there is no excitation of the diffusion mode. We also observe 
in both cases the sound wave, that is more prominent in the second case. Note also that in both 
cases we have positive and negative velocities, that is, matter moving outward and inward (as 
it should be due to conservation of matter). 

5. Spectrum 

We use the previous hydrodynamic fields to calculate the final spectrum induced by the jet. To 
do so, we use the standard Cooper-Fry prescription. As our initial medium is static, we use 
fixed time freeze-out (neglecting the effect of the perturbation on the freeze out): 

dN _ f ^r-T7tZ+ 5 r E 6T , & 



d *Ve , 5 = t^j— + ^— , (15) 
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Figure 2. Velocity field v x created by a jet moving along the x direction (c^ = 1/3, 
T s = l/(47rT), a = r s ). The jet is assumed to disappear at t = 7 fm while the spectrum 
calculated at t = 10 fm. The two figures (a) and (b) are for scenarios 1 and 2, respectively. The 
values of the parameters are arbitrary. Note that in (a) matter moves preferentially along the x 
direction, while in (b) it is in the Mach direction flljl. 



where V is the volume of the fireball and Tf rzt is the freeze-out temperature. According to our 
approximation, we expand the exponent to first order in the perturbation. We can now discuss 
two opposite regimes in the spectrum: 

Low energy particles E ~ Tf rzt . In this region we can expand the exponential in © and 
express the spectrum in terms of the energy and momentum deposited (to first order): 

dN -*.[ E E dep P P dep \ 

— = e /«* V + — v — + -f- . (16) 

d 6 p \ Tf rzt e + Po Tfrzt £o + Po J 

One finds that soft particles are insensitive to the particular shape of the flow field and their 
angular dependence is just a cosine of the relative angles of the observed particle and the jet. 

High energy particles E » Tf rz t which compensates the smallness of the flow velocity. 
The integral is now dominated by the maximum of the exponent. Thus, only the points of 
maximum modifications of hydrodynamic fields contribute to the final spectrum. It is clear 
then that in the scenario 1 the diffusion mode totally dominates the spectrum and there are no 
angular correlations related to the speed of sound. However, in the second scenario one does 
not excite that mode, and the final spectrum do reflect the shape of the sonic disturbance. 

In order to illustrate the effect of the modification of the final particle production due to a 
jet moving with y = and (j) = ir, we show in FigEJthe normalized spectrum defined as follows: 

1 dN r 2n dN 

c= Qi^^ =a ^ ) ■ Q =L d ^ (!/=(W> - (17) 

We observe that, as claimed, the effect of the diffusion mode hides the Mach cone in the final 
production. As a consequence, we do not observe such a structure in fig FigOJ^a); we do observe 
it, however, in FigOJb), where such mode is not excited. 

The parameters in FigEJare arbitrary (with the exception of T s , which is set to its minimal 
bound ^2])) and are chosen as a matter of illustration. Let us remark that these spectra are 
very sensitive to the value of T s , which could provide a experimental constraint on its value, 
provided a correct understanding of the deposition process. 



° 0.4 








0.38 




— Pt= T frzt 








....p t =10 T llzt 




0.36 




Pt=20 T M 




0.34 








0.32 








0.3 








0.28 












i i i I ■ i i i I i i i i I 





O 1.4 F 
1.2 - 
1 

0.8 
0.6 
0.4 - 
0.2 - 



— Pt=T, rzt 
....p t =10T fr 

P,=20T (r . 



Figure 3. The normalized spectrum of associated secondaries versus the azimuthal angle (p. 
Three curves are for different p t at y = for c 2 s = 1/3, T s = 1/(4ttT), a = T s . Note the different 
scales. The jet disappears completely at t = 7 fm while the spectrum is calculated at t = 10 fm. 
The two figures (a) and (b) are for scenarios 1 and 2, respectively. 



6. How can the effect be observed? 

The most important feature of conical flow is its direction, which is normal to the shock front 
and thus making the so called Mach angle (J2) with the jet direction, determined basically by 
the speed of sound in matter. So one can conclude that quenched jets must be accompanied by 
a cone of particles with the opening angle 9m- This cone angle should be the same for any jet 
energy (in contrast to radiation angles, which are shrinking at high energy). 

The appropriate value for the speed of sound is a time-weighted average of three stages: (i) 
the QGP phase (c s ~ l/y3), (ii) the mixed phase c s ~ and (hi) the hadronic or "resonance 
gas" stage, with c s as \A2 ^3]. For RHIC we found the time- weighted average to be c s as .33. 
We thus conclude that in this case the emission angle of the conical flow should be at angles 
about 70 degrees relative to the jet, or in radians 

demission = ^ ± arccos(c s /c) 1.9,4.3. (18) 

The first measurements of interest to this issue have been made by the STAR collaboration 
[Tlj . H5| which studied two particle correlations in which the trigger particle has 4 < p t < 6 GeV ', 
while the associated particles have basically all momenta (in fact 0.15 < pt < 4 GeV ). At 
Ac/) = one finds a peak due to particles from the triggered jet, while particles from the 
companion jet should result in a peak at A<j) = ir. Such peak is clearly seen in pp collisions, in 
which no matter is present: but in Au-Au collision one finds instead a double-peaked distribution 
with a minimum at the original jet direction A<j) = tt. Remarkably this value it ± 1.2 nicely 
agrees with the observed positions of the two maxima of the distribution of secondaries. 

After submission of our paper in the preprint form, two more important observations have 
been made, both discussed at this workshop [T5].jl6j. 

The PHENIX collaboration have shown JH] correlation function in which the companion 
particle is harder than the average pt > 1 GeV. They also reported data for 6 centrality 
classes and in all of them (but the most peripheral one) there is a minimum at 180° and even 
sharper peaks at the same angle as we predict. We also remark that our results in Fig. 3(b) for 
Pt = WTf rzt (Tf r zt = 100 MeV) are similar to the correlation functions in [TT)] . 

The STAR collaboration also reported ^H] the associated particle mean (pt) (</>), which also 
displays maxima away from <p = tt, where one finds a clear minimum. This fact strongly 
supports our explanation of the peaks by conical flow. The fact that at 180° the mean pt of 
associated secondaries is consistent in magnitude with the mean transverse momentum in the 



average background events means that in the experimental conditions the inward-moving jet is 
completely quenched by freeze-out and no contribution from the "nonhydro core" is visible in 
the data. 

Unfortunately, we can only see a projection on Acp because the distribution in rapidity of 
the associated jet is very wide, wiping out another projection of the cone. One could think that 
this wide rapidity distribution could also wipe out the double peak structure in the final two 
particle correlation. Even though the full hydrodynamic problem should be solved, we can give 
here an argument by which this does not happen. For that we assume boost invariance of the 
medium and that the rapidity (yj) distribution of the away jet is flat. We will also assume that 
particles are formed in a cone around the jet axis. It is clear then, that in a frame with rapidity 
yj (where the away side jet is transverse) the spectrum of those particles can only depend on 
the momentum of the particle (pj) and its angle with respect to the jet (9*). We assume also 
that this dependence factorizes and that the dependence in 9* is a very narrow function of angle 
9* (the cone opening angle). Thus, the spectrum of particles produced by a jet with rapidity yj 
can be written as, 

dN 

where P(yj) is the probability of finding a jet with such rapidity. We will concentrate now 
on the spectrum of particles at mid rapidity (y=0) in the lab frame. Boosting back to this 
frame, we find pj = pt cosh(j/j) and cos(</>) = cos(6>*) cosh(yj), where pt and 4> are the transverse 
momentum and azimuthal angle respectively. Thus: 



A-(cos(n - c-os(9 c *)) = — , COsh( y *) = . (20) 



S(yj - y 

cos 2 (<p) 

Upon integration on yj (assuming P constant) we find 

which is clearly a peaked distribution around <j) = ±6>*. What is more, f p should be a steeply 
falling function (exponential in our case) of its argument, which depletes more the fill up. A 
width on the angular dependence will modify the result but a quantitative answer to this effect 
requires further study. 



7. Discussion: 

We have considered an idealized case of homogeneous matter at rest. The determination of the 
exact shape of this cone for real collisions is not simple. It is of course just a technical matter 
to include the superposition of radial, elliptic and conical flows in a single hydro simulation. 
The open issue, however, would be the inclusion of viscosity at the hadronic stage, which is not 
supposed to be small. Presumably, as for elliptic flow, the use of a hadronic cascade afterburner 
would provide more realistic results. 

When the conical solution reaches the surface of the fireball or the freeze-out surface, the 
energy stored in outward (positive) flow would continue to move in the same direction, while 
the inward (negative) one would be absent. As the shock moves into the lower density region 
near the edge, energy conservation would force its amplitude to grow, like do the ocean waves 
approaching the beach, or the whip motion near its thin end. Furthermore, if the jet is energetic 



enough to punch through, the regions where cones reach the edge consist of two separated 
regions (circular rings if the jet goes via the diameter) moving toward each other and eventually 
colliding, where the cone and the fireball surface are tangent. 

We have seen that the choice of initial condition is very important for the final result, as the 
excitation of the diffusion modes may hide the cone formation in the final spectrum. We have 
presented two different scenarios that are, in a sense, opposite limits, as in one of then we do 
not excite this mode. The appearance of such mode in the actual initial conditions, will lead 
to some (or even complete) fill up of the cone in the two particle correlations. This fill up, is 
also expected for high energy particles, in which the jet punches through. This is however a 
completely different mechanism from the previous, which will lead to a very different j>t-spectrum 
of the associated particles (exponential in the first, power law in the second) 

Can there be an alternative explanation of the high-angle maxima observed? Another two 
different mechanisms have been proposed: 

Large angle gluon radiation It has been argued that QCD radiation calculated with 
Landau-Pomeranchuck-Migdal effect included, may have a maximum at large angles, similar 
to the Mach angle we discuss. Those calculations, however, used static (infinitely massive) 
scatterers. A simple kinematical calculation reveals that in the limit in which the energy of the 
jet is much larger than the energy of the target, any angular distribution of radiation in the center 
of mass frame leads to a forwardly peaked distribution in the rest frame of the target. In spite 
of this, even if one could argue that the scatterers were heavy, the same calculations presented 
in JJj cannot reproduce a double peak structure. In fact, the gluon distribution obtained in 
|17j are such that after integrating over the (flat) rapidity distribution of the associated jet, any 
double peak structure would disappear in the final correlation function. 

Deflection of the jet due to transverse expansion. This idea has been suggested in ^B] and 
is similar to those in ^H]- In this scenario, as the jet travels through the medium it changes 
its trajectory according to the direction of expansion of the medium it travels through. Thus, 
depending on the formation point of the jet, it will be deflected in a different direction with 
respect to the original axis (set by the triggered particle). Even though there are not quantitative 
estimates of this effect, we can establish one fundamental difference with respect to our scenario: 
particles are not produced in a cone around tt, as jets are deflected in either direction with respect 
to the original one. Thus, in an event by event basis, the two particle correlations should peak 
either at angles smaller or bigger than tt, and is the convolutions of all events which are the 
responsible for the double peak. This picture should be easily distinguished from ours via three 
particle correlations [20| . 

8. Summary: 

We suggest that the lost energy of the quenched jets is not just absorbed by the heat bath, but 
appears in the form of hydrodynamical collective motion similar to known "sonic booms" in the 
atmosphere behind the supersonic jets. The reason for that is that QGP seems to be a near- 
perfect liquid, with remarkably small dissipative effects and robust collective flows. We argue 
that this should result in a cone of particles moving in the Mach direction (^Q). The first data 
on the particle distribution associated with the quenched jet (the away side from high energy 
hadron trigger) indeed show two peaks, with cone angles which agree well with our prediction. 
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